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Abstract 

A power law is postulated for both site and bond percolation thresh- 
olds. The formula writes p c = po[(d — l)(q — l)]~ a d b , where d is the 
space dimension and q the coordination number. All thresholds up to 
d — > oo are found to belong to only three universality classes. For first 
two classes b = for site dilution while b = a for bond dilution. The last 
one associated to high dimensions is characterized by b = 2a — 1 for both 
sites and bonds. Classes are defined by a set of value for {po; a}. De- 
viations from available numerical estimates at d < 7 are within ±0.008 
and ±0.0004 for high dimensional hypercubic expansions at d > 8. The 
formula is found to be also valid for Ising critical temperatures. 
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Percolation theory has been known for several decades [1, 2]. It deals 
with the effects of random dilution of either sites or bonds in a lattice. Flory 
introduced it first in the framework of chemical industry [3]. Upon site or bond 
dilution a sharp change is found to occur in the connectivity of the system at 
some threshold p c in the density of occupied site or bond. Although it is a 
purely geometrical phenomenon, this change can be described in terms of usual 
second order phase transitions. This mapping to critical phenomena made 
percolation a full part of the theoretical framework of collective phenomena and 
statistical physics. It is indeed a very powerful and general tool. Percolation 
has thus been applied to numerous problems in a large variety of fields, even 
outside physics (for reviews, see [4, 5, 6, 7]). 

However, despite both its success and its mathematical ground, perco- 
lation theory has resisted exact calculations. Most known data are numerical 
estimates, from both Monte Carlo simulations and series expansions [6, 8]. In 
particular, analytic calculations of percolation thresholds has proven to be a 
rather difficult task. For instance, twenty years or so were necessary to prove 
the numerical estimate of p c = 1/2 for the square bond percolation threshold 
[9]. The bond threshold is also known exactly for the two-dimensional hon- 
eycomb and triangular lattices [9]. The situation is even worse in the case 
of site percolation. Indeed, thresholds are known exactly only in the cases of 
two-dimensional triangular and Kagome lattices [9]. At dimension higher than 
two, no thresholds were determined exactly. From simulations, percolation 
thresholds are found to depend on both the space dimension d and the coor- 
dination number q. Only the pathological Cayley tree has been solved exactly 
to yield p c = l/(q — 1) [10]. This expression which holds for both bonds and 
sites, yields good results at high dimensions. 

Along the general form of the Bethe expression, various empirical for- 
mulas have been tried to yield large classes of thresholds without success. For 
instance the formula p c = d/[(d — l)(q — 1)] was proposed twelve years ago for 
site percolation [11]. Another one, p c = d/[(d — l)q] was suggested for bond 
dilution [12]. Unfortunately, both formulas are not satisfactory for all lattices, 
especially at low dimensions. More recently the expression p c = l/y/q — 1 
was derived for site percolation thresholds [13, 14]. Very good results are 
obtained at two dimensions but not at higher dimensions. Another similar 
form was also noticed elsewhere [15]. In this paper, for the first time, one 
unique power law is found to yield within an excellent accuracy both site and 
bond percolation thresholds for all lattices at all dimensions. The power law 
is p c = po{{d — l)(q — l)}~ a d b . From a log-log plot all available data (d < 7) 
are found to fit on two straight lines. In both cases b = for site dilution and 
b = a for bond dilution. One line includes two-dimensional triangle, square 
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and honeycomb lattices, which constitute the first class, characterized by 
{po = 0.8889; a = 0.3601} for site dilution and by {p = 0.6558; a = 0.6897} 
for bond dilution. Two-dimensional Kagome and all other regular lattices (for 
d > 3) align on the other unique line and constitute the second class, char- 
acterized by {po = 1.2868; a = 0.6160}, and {p = 0.7541; a = 0.9346} 
for sites and bonds respectively. Deviations from available numerical esti- 
mates are very small and range from to ±0.008. At high dimensions a third 
class is found which recovers the infinite Cayley tree limit. It is defined by 
{b = 2a — 1; p = 2 a ~ 1 } for both sites and bonds. Using high dimensional 
hypercubic expansions [14, 17, 18] as data we found {a = 0.8800} for sites and 
{a = 0.3685} for bonds. The agreement is excellent with deviations within 
±0.0004. The value of dimension d c at which the third class holds is between 
d — 6 and d — 8. The results and perspectives are discussed. 

From an analysis of the Bethe approximation on a regular lattice, we 
found recently [13, 14] an underlying lattice homogeneity breaking. On this 
basis we obtained a new expression, 

x c = l/yj(q - 1) , (1) 

for the percolation threshold, instead of the well known Bethe expression, 

*? = V(9-1). (2) 

Eq. (1) yields results which are in good agreement with exact and nu- 
merical site percolation thresholds at d — 2. However results get poor at higher 
dimensions. This discrepancy evidences a missing d-dependence which could 
be equal to one at d — 2. It thus hints a possible (d — 1) dependance. In 
parallel Eq. (1) did single out a (q — 1) dependence. From the hypercube case 
we could expect dimension and coordinance to play a similar role. Therefore 
one extension of Eq. (1) is to substitute the product (d — l)(q — 1) for (q — 1). 
On this basis, a natural generalization of Eq. (1) is the power law, 

p c = p [(d-l)(q-l)]- a . (3) 

It is worth to note that Eq. (3) reduces to Eq. (1) for the set {po = 1; a = 1/2} 
besides the variable multiplicator factor (d — 1) which equals one at d — 2. 
At this stage, we report a log-log plot of known estimates p c [7, 19] as a 
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function of the associated product (d — l)(q — l). The plot is shown in Fig. (1) 
for the site case. All the points align on only two straight lines, corresponding 
to the two classes of lattices, above mentioned. 

In parallel the log- log plot for bond data [7, 19] exhibits constant de- 
viations from a straight line at each dimension. These deviations suggest a 
dimension rescaling of the variable (d — l)(q — 1) resulting in, 

The log-log plot using Eq. (4) is shown in Fig. (2). The situation is now 
identical to the site case, as all points align on two straight lines associated 
to the two classes of lattices. The percolation thresholds using Eqs. 3 and 
4 are shown in Tables (1) and (2) together with exact results and numerical 
estimates for respectively site and bond dilution, respectively. 

At this stage we have a one exponent power law for sites (Eq. 3) and 
bonds (Eq. 4). However the site variable (d — l)(q — 1) has to be rescaled by 
dimension in the case of bonds. The two expressions can thus be unified under 
the form 

p c = Po [(d-l)(q-l)]- a d b . (5) 



From Eqs (3) and (4) we have b = for sites and b = a for bonds. The 
agreement is remarkable and holds true from d = 2 up to d = 7 where threshold 
estimates are available. 

To complete the discussion we compare our results to d-dimensional sim- 
ple hypercubic lattice percolation thresholds derived as l/(q — 1) expansions. 
For site percolation it is [16], 

P? = (Q ~ + |(? - I)" 2 + - I)" 3 + ™(Q- I)' 4 + , (6) 
which becomes for bond percolation [17], 

pf = (<Z - I)" 1 + \{q~ I)" 3 + y (? - I)" 4 + 57(g - l)" 5 + ... , (7) 
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In parallel from Eq. (5) we get the expansion, 



p c = p 2«- b ( q -l) b -^{l + (a+b)(q-l)- l + - + a + 2 f + * \ q-ir + ■■■}■ 

(8) 

Eqs. (6) and (7) have the same first leading term l/(q— 1) which is the 
Bethe result. In our case we have, 

p 2 a -\q-l) b - 2a , (9) 

which becomes p c = po2 a (q — l) _2a for sites with the exponent 2a = 1.23 and 
Pc = Po(q — l) -a for bonds with a = 0.94. Both cases are clearly different 
from the Bethe result l/(q — 1). The respective l/(q — 1) exponent is close 
to, but definitively different from one. Therefore we do not recover the Bethe 
asymptotic limit. Moreover the site-p c becomes smaller that the bond-p c at 
d — 11 which is clearly non-physical. 

On this basis, although our formula agrees perfectly with available nu- 
merical data up to d = 7 (see Tables 1 and 2), a third class must exist at high 
dimensions. To determine its characteristics we first note that the Cayley tree 
result — l/(q — 1) is believed to be the exact d — > oo asymptotic limit for 
both bonds and sites. It then indicates that two different constraints 6 = 
and b = a for respectively sites and bonds cannot hold at high dimensions. 
The 1/(5 — 1) limit is recovered from Eq. (9) if, and only if b = 2a — 1 together 
with p = 2 a ~ l . Eq. (5) becomes, 

p c = 2 a - 1 [(d-l)(q-l)}~ a d 2a -\ (10) 

which gives a straight line in a log-log plot of 2dp c versus 2d 2 /[(d — l)(q — 1)]. 

We then determine the value of exponent a using numerical estimates 
from the l/(q — 1) expansions in Eqs. (6) and (7) which are supposed to be 
exact in the d — > 00 limit. From Fig. (3) we find a = 0.8800 for sites and 
a = 0.3685 for bonds. It is worth noting that while Eqs. (6) and (7) arc 
derived for hypercubes only, Eq. (10) holds for any lattice at d > d c . 

From Fig. (3) the best estimate for the crossover dimension is d c — 8. 
However we cannot preclude that d c is as small as d c = 6 which is the upper 
critical dimension for percolation. Differences between numerical estimates 
and associated values of p c deduced, on the one hand from Eq. (3) or (4), 
and on the other hand from Eq. (10), are comparable to numerical errors 
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(see Table (3) and Fig. (3)). From current data we have 6 < d c < 8. Exact 
determination of d c requires more accurate numerical estimates of percolation 
thresholds at these dimensions. 

Moreover the crossover at d c seems to be driven by a dimensional phe- 
nomena. This conjecture is based on the following observation. At d — 5, fee 
lattice which has (d — l)(q — 1) = 156 belongs to the second class (see Fig. 
(1)). In parallel hypercube at d = 9 belongs to the third class (see Fig. (3)) 
though it has (d — l)(q — 1) = 136. 

The various universality classes for site and bond dilution were found 
to be identical. One class is restricted to two dimensions and contains trian- 
gle, square and honeycomb lattices. The second universality class embodies 
Kagome and all lattices at 3 < d < d c . It is interesting to stress that square 
and Kagome which both have d = 2 and q = 4 turn out to be in different 
classes. The third class embodies all lattices at d > d c . At d — 1 (q — 2) exact 
value p c = 1 for both sites and bonds implies a = 0. 

At this stage it is worth to emphazise numerical values of po and a have 
been determined using as input data what is denoted in Tables (1) and (2) as 
"exact" thresholds. Obviously these values will be different using a different 
set of input or a restricted one. We chose arbitrarily to use and report values 
form Refs. (6, 17). We are here advocating the power law form of Eq. (5) 
rather than the third decimal associated to the determination of po and a. 

Error bars for numerical estimates are believed to be on the last given 
figure. It should thus not exceed ±0.001 in case p c is given with four decimals 
and ±0.01 in case it is given with three decimals. Therefore deviations A = 
Pc — Pc are significant at d = 2 though they are small. For all lattices at this 
low dimension, | A | exceeds 0.001 (A is only 0.0013 in the triangular lattice 
case, but there, p c is known exactly). 

To further investigate this anomaly at d — 2, we have explored percola- 
tion thresholds when the connection range includes next (nnn) and next-next 
(nnnn) neighbors. Numerical estimates pi are scarce and restricted to site 
percolation [20]. They strongly suggest all d = 2-lattices enter the second 
universality class as soon as the connection range exceeds nearest neighbors. 
Associated site percolation thresholds are then calculated using for q the sum 
of sites within the connection range. We get respectively p c = 0.388 for nnn- 
square, 0.294 for nnn-triangular, nnnn-square and honeycomb, and 0.225 for 
nnnn triangle. It means 10 3 A = —1, ±1, +2, —6, respectively, which is 
always smaller than the numerical estimate error bar ±0.01. However available 
square lattice site percolation thresholds for longer ranges [21] don't belongs 
to this class. It shows that an additional class should be added to account for 
long range interaction percolation (over nnnn). 
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At d > 3, errors A may not have a statistical significance. Only the 
value of the bond percolation threshold for the simple cubic lattice exceeds 
significantly 0.001. Moreover, the random distribution of deviations opposes 
the existence of a missing systematic correction to Eq. (5). 

Last but not least we found that our power law is also valid to predict 
Ising critical temperatures T c . Indeed the phase transition nature of percola- 
tion makes critical thresholds similar to critical temperatures. Dealing with 
pair exchanges, it is natural to use the bond percolation formula (Eq. 4) for 
the reduced temperature K c = j^r, with J the exchange coupling. Numerical 
estimates for K c at d = 3 are K e c = 0.2217, 0.1575, and 0.1021 for sc, bcc 
and fee lattices respectively [22]. The set {p = 0.6525; a = 0.9251} in Eq. 
(4) leads to critical temperatures K c which depart from these values by the 
amount K c - K% = +0.0012, +0.0008, and -0.0006, respectively. It is again 
comparable to the numerical estimate errors. 
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Figure captions 



Figure 1: Inverse of site percolation thresholds as a function of 
(q — l)(d — 1) in logarithmic scales. 

Figure 2: Inverse of bond percolation thresholds as a function of 
(q — l)(d — l)/d in logarithmic scales. 

Figure 3: Percolation thresholds in high-d hypercubes using both numer- 
ical estimates for d = 6 and d = 7 and hypercubic expansions from Refs. (14) 
and (15) at d — 8, 9,..., 15, d = 20, 30, 50. The origin corresponds to the 
Cayley tree limit at infinite d. Straight lines are according to Eq. (10). Broken 
curves correspond to Eqs. (3) and (4), for the second universality class. 
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Table 1: Site percolation thresholds from this work p c compared to "exact esti- 
mates" p e c taken from [6, 17]. A = p c — p e c . * means not included to determine 
Po and a. The first universality class is defined by p = 0.8889 and a = 0.3601. 
The second universality class is defined by p = 1.2868 and a = 0.6160. 
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Table 2: Bond percolation thresholds from this work p c compared to "exact esti- 
mates" p e c taken from [6, 17]. A = p c — p e c . * means not included to determine 
p and a. The first universality class is defined by p = 0.6558 and a = 0.6897. 
The second universality class is defined by p = 0.7541 and a = 0.9346. 
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Table 3: High-d hypercube thresholds: p c -site and p c -bond from Eq. (10) and the 
"hypercube" expansions p s c and p\ from respectively [14] and [15]. At d > 10 the 
difference is smaller than 10~ 4 ; even relative errors are negligible, as it can be seen 
in Fig. 3. 
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